# Aufgabe 7.1

#Konstanten vereinbaren
a <- 2
e <- 1/a

# Stichprobenumfang
n <- 100

# Anzahl der Stichproben
s <- 1000

# Berechnung der Schätzungen
y <- c(rep(s,0))
for(i in 1:s)
{x <- runif(n)^e;h<-mean(x) ;h <- 1/h-1; y[i] <- 1/h}

# Auswertung
ss <- summary(y); ss

par(mfrow=c(2,1))

# Histogram
hist(y, main="Verteilung der Momentenschätzung für den 
              Parameter a=2 bei s=1000 Stichproben")
# Boxplot
boxplot(y, main="Boxplot der Momentenschätzung")